Transmission and Reflection of Bose-Einstein Condensates Incident on a Gaussian 

Potential Barrier 



A.M. Martin 1 , R.G. Scott 2 ' 3 and T.M. Fromhold 3 
1 School of Physics, University of Melbourne, Parkville, Victoria 3010, Australia 
2 Department of Physics, University of Otago, P. 0. Box 56, Dunedin, New Zealand 
School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, United Kingdom 

(Dated: February 6, 2008) 

We investigate how Bose-Einstein condensates, whose initial state is either irrotational or contains 
a single vortex, scatter off a one-dimensional Gaussian potential barrier. We find that for low atom 
densities the vortex structure within the condensate is maintained during scattering, whereas at 
medium and high densities, multiple additional vortices can be created by the scattering process, 
resulting in complex dynamics and disruption of the atom cloud. This disruption originates from 
two different mechanisms associated respectively with the initial rotation of the atom cloud and 
the interference between the incident and reflected matter waves. We investigate how the reflection 
probability depends on the vorticity of the initial state and on the incident velocity of the Bose- 
Einstein condensate. To interpret our results, we derive a general analytical expression for the 
reflection coefficient of a rotating Bose-Einstein condensate that scatters off a spatially-varying 
one- dimensional potential. 
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I. INTRODUCTION 

Recently, several experiments have investigated the 
scattering of Bose-Einstein condensates (BECs). These 
experiments have included Bragg reflection in an opti- 
cal lattice , reflection from optical 2] and magnetic Q 
mirrors, diffraction from a grating Q and quantum reflec- 
tion from a silicon surface In each case, interest has 
focused on the reflected component of the BEC. For ex- 
ample, investigations of quantum reflection from a silicon 
surface have revealed that inter-atomic interactions have 
a dramatic effect on the internal structure of the atom 
cloud [5|, |6|. So far, reflection experiments have been 
restricted to condensates whose initial state contains no 
dynamical excitations. However, the methodology for 
creating and observing vortices in BECs is well estab- 
lished 

EE II E3 and, in the case of Bragg reflection, 
numerical simulations predict that vortices and solitons 
in the BECs initial state strongly influence the subse- 
quent dynamics Previous theoretical work has also 
shown that the presence of a vortex in the initial state 
can have a pronounced effect on the internal structure 
of a BEC that undergoes classical reflection from a hard 
wall atom mirror [l2,[l3[. To our knowledge, however, 
there has been no consideration of the effect of vortices 
on BECs approaching a potential barrier of finite width 
and height, where quantum- mechanical tunneling is pos- 
sible as well as reflection. Quantum tunneling of BECs 
can be studied experimentally by using sheets of laser 
light to create the potential barrier and plays a cru- 
cial role in the dynamics and macroscopic coherence of 
cold atoms in optical lattices [l5(- 

In this paper, we investigate how BECs scatter off a 
Gaussian potential barrier, which allows both reflection 
and quantum- mechanical tunneling of the atoms. We use 
numerical simulations of the Gross-Pitaevskii equation to 



make a detailed study of how the strength of the inter- 
atomic interactions, the energy of the incident BEC, and 
the vorticity of the initial state affect the scattering prop- 
erties of the condensate. Our simulations reveal regimes 
in which dynamical excitations disrupt both the reflected 
and transmitted atom clouds. In one regime, which we 
call rotational disruption, the excitations originate from 
the effect of the initial vortex on the scattering dynam- 
ics. By contrast, the regime of interferential disruption 
occurs both in the presence and absence of a vor- 
tex in the initial state with the excitations being cre- 
ated from interference between the incident and reflected 
matter waves. We find that rotational disruption arises 
when the time taken for the BEC to scatter is compara- 
ble with, or exceeds, the rotational period of the vortex. 
Interferential disruption occurs when the scattering time 
is greater than the correlation time of the BEC, which is 
a measure of how quickly the atom cloud responds to a 
perturbation. To interpret our results, we derive a gen- 
eral expression for the reflection probability of a BEC 
containing a vortex in terms of the reflection probability 
of an irrotational BEC. We find that the vortex changes 
the reflection probability by altering the distribution of 
incident velocities at the barrier, due to the increase in 
the physical size of the cloud and, more importantly, the 
circulation of the atoms. 

To investigate the regimes of rotational and interfer- 
ential disruption we consider three BECs whose initial 
state contains a single vortex: BEC A with low atom 
density, BEC B with medium density and BEC C with 
high density. We identify the effect of the initial vortex 
by comparing the dynamics of these atom clouds with 
irrotational counterparts labeled BECs A{, Bi and C{. 

The layout of the paper is as follows: in section II we 
specify the parameters for each of the three BECs and 
describe the theoretical model and computational tech- 
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niques. In section III, we present our numerical results, 
which show how low, medium, and high density BECs 
scatter off the Gaussian potential barrier, both with and 
without a vortex in the initial state. In section IV, we 
derive a general expression for the reflection probability 
of an irrotational BEC. We then use this expression to 
derive an approximate analytical formula for the reflec- 
tion probability of a condensate whose initial state does 
contain a vortex, impinging on the same scattering po- 
tential. We compare this formula with reflection proba- 
bilities obtained directly from the numerical simulations 
presented in section III and use it to provide physical 
insight into the effect of vortices on the scattering pro- 
cess. Finally, in section V, we summarize our results and 
propose experiments to test them. 



II. SYSTEM PARAMETERS, THEORETICAL 
MODEL, AND COMPUTATIONAL TECHNIQUES 

Each BEC contains N 23 Na atoms of mass m and 
is initially confined by a harmonic trapping potential, 
Vt(x, y, z) = m[uo 2 (x + Ax) 2 + u 2 y 2 + co z z 2 ]/2 centered 
at (—Ax, 0,0). We consider trap frequencies uo x = uo y = 
50 x 2tt rad s _1 < w z so that the spatial width of the 
BEC is much smaller along the z-direction than along the 
x- and ^/-directions. Consequently, the dynamics reduce 
to two-dimensional motion in the x — y plane. At time 
t = 0, we create an additional Gaussian potential barrier 



V L (x) = V e~ 



(1) 



of width <j = l /im along the x-direction, by switching 
on a far blue-detuned laser beam that travels along the 
^/-direction and creates a sheet of laser light. The in- 
tensity of the laser beam determines the barrier height, 
Vo, which we take to be 6.2 peV, similar to recent ex- 
periments Simultaneously, we accelerate the BEC 
towards the Gaussian potential by abruptly displacing 
the harmonic trap through a distance Ax along the re- 
direction 0, H, @- After the displacement, the center 
of the trap coincides with the Gaussian potential energy 
barrier and the total potential energy of the trap and 
laser beam, for motion in the x — y plane, is given by 



V(x,y) 



m(u 2 x x 2 + u 2 y y 2 )/2 + V L (x). 



(2) 



The solid curve in Fig 1(a) shows the form of V(x, y = 0). 
After the displacement of the harmonic trap, the conden- 
sate moves away from its initial state [shown schemati- 
cally by the dashed curve in Fig. 1(a)], reaches the Gaus- 
sian potential with a mean incident velocity v x w u x Ax, 
and is then scattered by the potential barrier. 

The time-dependent Gross-Pitaevskii equation for the 
system is 



ih 



dt 



'2m 



dx 2 



dy 2 
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FIG. 1: (a) Solid curve: A schematic of the potential energy 
V (x,y = 0) created by the harmonic trap and laser beam 
when t > 0; Dashed curve: a schematic of the initial proba- 
bility density |^(x, y = 0, t = 0) | 2 for BEC B. (b) The initial 
density \^f(x,y,t = 0)| 2 (black=high density, white=zero) of 
BEC B, where the white arrow shows the direction of the 
condensate's circulation and the horizontal bar denotes the 
scale, (c) The equivalent phase plot, <fi(x, y, t = 0), (white= 0, 
black= 2tt). 



[V(x,y) + U 



9(x,y,t) \ 2 ]^(x,y,t) 
(3) 



where Sfr (x,y,t) is the wavefunction for motion in the 
x — y plane at time t > and Uo = 47rh 2 a/m where the 
s- wave scattering length a = 2.9 nm. 

We determine the initial BEC wave function by solv- 
ing Eq. (3) for t < using an imaginary time algorithm 
[l7j . When the initial state contains a vortex, we im- 
pose the requirement that there is a 27r-phase change 
in the condensate wavefunction around the trap center 
at (x,y) = (— Ax,0), which corresponds to a quantized 
angular momentum of — % about the z-axis. The wave- 
function is normalized according to 



/ 



^ (x,y,t) | 2 dxdy = 



(4) 



where L z is the confinement length in the z-direction. 
For BECs A, B and C we choose N/L z as 2.5 x 10 8 
m _1 , 5 x 10 9 m _1 and 2.5 x 10 11 m _1 respectively, which 
gives corresponding peak atom densities no = 2.8 x 10 18 
m -3 , 2.1 x 10 19 m -3 and 1.6 x 10 20 m -3 . The density, 
| ty(x,y,t = 0) | 2 , and phase, <fi(x,y : t = 0), of the ini- 
tial state for BEC B are shown in Figs. l(b,c). Having 
obtained the initial state of the BEC, we determine its 
motion by solving the Gross-Pitaevskii equation (3) nu- 
merically using the Crank-Nicolson method |l8j . 
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FIG. 2: Evolution of BEC A: plots of |^(x, y, t)\ 2 (black=high 
density, white=zero) for v x — 6.3 mm s _1 at t = 3.6 ms (a), 
4.0 ms (b), 4.4 ms (c), 4.8 ms (d), 5.2 ms (e), 5.6 ms (f), 6.0 ms 
(g) and 6.4 ms (h). Coordinate axes are inset and the hori- 
zontal bar indicates the scale. The arrows in (a) and (h) show 
the direction of rotation. 

III. NUMERICAL RESULTS FOR SCATTERING 
OFF A GAUSSIAN POTENTIAL BARRIER 



In Fig. 2, we show the density profile of BEC A, with 
Ax — 20 fim (v x = 6.3 mm s -1 ) at t = 3.6ms (a), 4.0ms 
(b), 4.4ms (c), 4.8ms (d), 5.2ms (e), 5.6ms (f), 6.0ms 
(g) and 6.4ms (h). Figs. 2(a-e) show that as the BEC 
approaches the Gaussian scattering potential at x = 
(dashed lines in Fig. 2), a standing wave forms due to 
interference between the incoming and reflected matter- 
waves. In Figs. 2(d,e) the standing wave undergoes a tt 
phase shift between the upper (large y) and lower (small 
y) edges of the BEC. This is due to the no n- uniform ini- 
tial phase of the BEC, shown in Fig. 1(c). After scatter- 
ing, the BEC splits into reflected (x < 0) and transmitted 
(x > 0) components [Figs. 2(f-h)]. Each of the compo- 
nents contains a vortex [enclosed by the arrows in Fig. 
2(h)]. For the transmitted cloud, the quantized circula- 
tion is in the same direction as in the incident cloud. By 
contrast, the reflected cloud rotates in the opposite direc- 
tion 

HD3- Physically, this is because atoms towards 
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FIG. 3: BEC A: plots of \V(x,y,t = 10ms)| 2 for v x = 
9.4 mm s -1 (a), 7.5 mm s -1 (b) and 6.3 mm s -1 (c). The circu- 
lar arrows indicate the direction of rotation and the horizontal 
arrows show the values of Lr and Lt, defined in the text and 
footnote 19]. Co-ordinate axes are inset. 

the top of the rotating incident cloud, in Fig. 2, approach 
the barrier with a higher velocity component along the 
x-direction (v x ) than those at the bottom. But because 
v x is reversed after reflection, the direction of rotation is 
also reversed in the reflected cloud. This is in contrast to 
the transmitted atoms, which emerge through the barrier 
with little change in their velocity and so the direction 
of rotation is preserved. Since the trap remains switched 
on throughout the scattering process, we do not observe 
the splitting of the reflected component of the BEC that 
was reported in Ref. |l3j . 

We now consider how changing v x affects the dynam- 
ics of the atom cloud. Figure 3 shows the density profile 
of BEC A at t = tt/uj x = 10 ms after trap displacements 
of 30 jam (a), 24 jam (b) and 20 jam (c), corresponding 
to v x = 9.4 mms -1 , 7.5 mm s _1 and 6.3 mm s -1 respec- 
tively. Since this time is half the period of oscillation 
of the trap, both the transmitted and reflected portions 
of the condensate will approximately be at their turning 
points, and hence nearly stationary. In Fig. 3, the ver- 
tical dashed line marks the trap center at x = 0. When 
v x = 9.4mms _1 , the average kinetic energy of the atoms 
incident upon the barrier is lOpeV, which exceeds the 
Gaussian barrier height. Hence, most of the condensate is 
transmitted [Fig. 3(a)]. As we decrease v x [Figs. 3(b,c)], 



(a) 




(b) 



(d) 



(e) 



y 



(c) 



-►X 



20|um 




FIG. 4: Evolution of BEC £: plots of |^(x, y, t)\ 2 (black=high 
density, white=zero) for v x = 6.3 mm s _1 at t = 3 ms (a), 4 ms 
(b), 5 ms (c), 6 ms (d), 7 ms (e), 8 ms (f). The dashed line at 
x — marks the point where the laser potential is maximal. 
Co-ordinate axes are inset and the horizontal bar indicates 
scale. Lower plot: phase <p(x,y,t = 8 ms) [white= 0, black= 
2tt] within the region enclosed by the box in (f). Arrows 
indicate the direction of circulation. 

the average energy of the incident atoms decreases and 
more of the condensate is reflected. For v x = 6.3mms _1 
[Fig. 3(c)], most of the atoms are reflected by the barrier 
because the average kinetic energy of the incident atoms 
is 4.6 peV, which is less than the barrier height. Since 
the BEC has a finite spatial width along the x-direction 
of l x = 13 //m, atoms that are towards the left-hand side 
of the BEC's initial state travel further before reaching 
the barrier and therefore have a higher incident veloc- 
ity. Such atoms have a higher transmission probability 
and therefore form a large fraction of the transmitted 
atom cloud. We therefore expect that the distance (Xt) 
that the higher velocity transmitted cloud travels past 
the scattering potential before coming to rest in the har- 
monic trap will be greater than the distance (Lr) that 
the slower reflected cloud retreats from the scatterer be- 
fore reaching the turning point of the harmonic trap . 
Our numerical simulations confirm this: for example, in 
Fig. 3, Lt > Lr for all the values of v x . 

We now consider how the higher density BEC B re- 
flects off the potential barrier. Figure 4 shows the den- 
sity profile of BEC B at t = 3 ms (a), 4 ms (b), 5 ms (c), 
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FIG. 5: Evolution of BEC Be plots of \^{x,y,t)\ 2 
(black=high density, white=zero) for v x — 6.3mms _1 at t = 
3ms (a), 4ms (b), 5ms (c), 6ms (d), 7ms (e), 8ms (f). The 
dashed line at x = marks the point where the laser potential 
is maximal. Co-ordinate axes are inset and the horizontal bar 
indicates scale. 

6 ms (d), 7 ms (e) and 8 ms (f), after a trap displacement 
of 20 /im (v x = 6.3mms _1 ). As the BEC impinges upon 
the Gaussian potential barrier, a standing wave forms be- 
tween the incoming and reflected matter waves. Figure 
4(b) shows the first stage of the standing wave forma- 
tion in which maxima (black) and nodal lines (white) 
appear at the leading edge of the atom cloud. In con- 
trast to BEC A, the reflected component of the BEC is 
significantly disrupted, as shown in Figs. 4(e,f). This 
disruption is accompanied by the formation of new vor- 
tices within the boxed region in Fig. 4(f). The phase 
variation, <fi(x,y,t = 8 ms), within this region is shown 
in the lower part of Fig. 4(f). Around each vortex, <\> in- 
creases from to 2tt in the direction of circulation shown 
by the arrows in Fig. 4(f). 

To explain the disruption shown in Fig. 4, we first con- 
sider the dynamics of BEC Bi. Figure 5 shows density 
profiles for this initially irrotational condensate at t = 
3ms (a), 4ms (b), 5ms (c), 6ms (d), 7ms (e) and 8ms 
(f) after a trap displacement 20 /im (v x = 6.3 mms -1 ). 
No dynamical excitations are produced in the transmit- 
ted or reflected clouds. We therefore conclude that the 
disruption observed for BEC B (Fig. 4) is related to the 
rotation of the cloud. 

To gain further insights into this disruption, we now 
examine how changing v x affects the transmission and re- 
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BEC B: 
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plots of \^f(x,y,t = 10ms)| 2 for v x — 
9.4mms _1 (a), v x — 7.5mms _1 (b) and v x — 6.3mms _1 
(c). The dashed line at x — marks the point where the laser 
potential is maximal. The circular arrows indicate the direc- 
tion of rotation and the horizontal arrows show the values of 
Lr and Lt, defined in the text and footnote [l^. Co-ordinate 
axes are inset. 



flection of BEC B. Figure 6 shows the condensate density 
at t = 10 ms after trap displacements of Ax = 30 /im (a), 
24 /im (b), and 20 /im (c), corresponding to mean incident 
velocities ofv x = 9.4 mm s -1 , 7.5 mm s -1 and6.3mms _1 
respectively. For v x = 9.4mms _1 , the behavior of BEC 
B [Fig. 6(a)] is qualitatively the same as BEC A [Fig. 
3(a)]: the transmitted cloud retains the vortex structure 
of the incident cloud and only a small fraction of atoms 
are reflected. As v x decreases [Figs. 6(b,c)], the reflec- 
tion probability rises, but the structure of the reflected 
cloud becomes increasingly disrupted. This disruption 
contrasts with BEC A, where for all v x values the re- 
flected cloud retains a well defined vortex structure (Fig. 
3). For BEC £>, we find that Lt > Lr (see values in Fig. 
6) just as for BEC A. However, for v x = 7.5 mm s -1 and 
v x = 6.3 mm s -1 , Lr is lower for BEC B [Fig. 6(b,c)] 
than for BEC A [Fig. 3(b,c)] due to the fragmentation 
of the atom cloud that occurs for BEC B only. This 
fragmentation transfers kinetic energy from the center-of- 
mass motion of the reflected cloud into internal vorticity, 
thus damping the oscillation. 



The effect of the vortex on the scattering process de- 
pends on the scattering time of the vortex core (t sv ) rel- 
ative to its rotation time (t r ). Since the diameter of 
the vortex core is approximately the healing length £ = 

7r^ 2 m/h. 



1 / V87rnoa, it follows that t, 
If the ratio 



2£/v x and t r 



m , . Ax 
—7i£Axuj x OC — = 
2h vno 



(5) 



is ^> 1, so that there is insufficient time for the vortex to 
rotate during the scattering process, the vortex will have 
little effect on the dynamics of the transmitted and re- 
flected atom clouds. Physically, this is because all parts 
of the BEC are incident on the Gaussian potential with 
a variation in the incident velocity, due to the rotation, 
that is small compared to v x . Conversely, if t r /t sv ^ 1 we 
expect that the rotation of the BEC will disrupt the re- 
flected component of the BEC. We refer to such a regime 
as rotational disruption. For the lower-density BEC A, 
tr/tsv ^> 1 for all three of the v x values considered in 
Fig. 3, and so the reflected atom cloud is not signifi- 
cantly fragmented. By contrast, for the higher density 
BEC B t r /t sv < 1, and hence rotational disruption oc- 
curs in the reflected atom cloud. Note that reducing v x 
from 9.4mms _1 [Fig. 6(a)] to 6.3mms _1 [Fig. 6(c)] 
causes t r /t sv to decrease, and thus increases the frag- 
mentation of the reflected cloud. 

Disruption of the reflected cloud is even more pro- 
nounced for the higher density BEC C. Figure 7 shows 
the density of this BEC at t = 1 ms (a), 2 ms (b), 3 ms (c), 
4 ms (d), 5 ms (e), 6 ms (f), 7 ms (g) and 8 ms (h) after 
a trap displacement of 30 /im (v x — 9.4mms _1 ) |2£|. As 
the condensate impinges upon the scattering potential, a 
standing wave forms between the incident and reflected 
matter waves. Figs. 7(b-d). This standing wave seeds 
solitons [y, |2Jl |23 •> which decay via the snake instabil- 
ity |2i| into vortex-antivortex pairs [Figs. 7(e,f)], thus 
strongly disrupting the internal structure of the cloud 
[Figs. 7(g,h)]. Consequently, when incident atoms sub- 
sequently pass through the barrier they produce some ir- 
regularity in the transmitted atom cloud. Note, that this 
irregularity is less pronounced towards the right hand 
edge of the transmitted cloud, which contains the atoms 
that passed through the barrier before solitons and vor- 
tices formed at negative x (x = is marked by the dotted 
line in Fig. 7). We emphasize that although BEC C is 
in a regime where rotational disruption occurs, the ini- 
tial vortex does not (in contrast to BEC B) cause the 
severe fragmentation shown in Fig. 7. To demonstrate 
this, Fig. 8 shows BEC d (identical to BEC C, except 
irrotational) scattering off a Gaussian potential at t = 
lms (a), 2ms (b), 3ms (c), 4ms (d), 5ms (e), 6ms (f), 
7 ms (g) and 8 ms (h) after a trap displacement of 30 jam 
(v x = 9.4mms _1 ). Comparison of Figs. 7 and 8 reveals 
that the dynamics are qualitatively the same irrespective 
of whether or not the BEC contains an initial vortex. 

To understand these results, we recall previous work 
ElEl El El] on the Bragg reflection of a BEC in an op- 
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FIG. 7: Evolution of BEC C: plots of |tf(x, y, t)\ 2 (black=high 
density, white=zero) for v x — 9.4 mm s -1 at t = 1ms (a), 
2 ms (b), 3 ms (c), 4 ms (d), 5 ms (e), 6 ms (f), 7 ms (g) and 
8 ms (h). Co-ordinate axes are inset and the horizontal bar 
indicates the scale. 



tical lattice. In Refs. [21|,[22| it was shown that at Bragg 
reflection, fragmentation can arise from the density and 
phase imprinting that accompanies standing wave forma- 
tion. When the correlation time, 



rn 



noa 



(6) 



is much less than the Bloch period (£#), this imprint- 
ing leads to the formation of solitons and vortices, which 
disrupt the atom cloud. For BECs C and C^, a similar 
disruption occurs when t c < where t s = l x /v x is the 
approximate duration of the reflection process. This ef- 
fect is described as interferential disruption [l6|, since it 
originates from the interference pattern, in this case pro- 
duced by the superposition of the incident and reflected 
matter waves. For BEC C this interferential disruption 
dominates the dynamics and completely masks any ef- 
fects due to the rotational disruption. 
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FIG. 8: Evolution of BEC Cn plots of \^(x,y,t)\ 2 
(black=high density, white=zero) for v x — 9.4 mm s -1 at t = 
lms (a), 2ms (b), 3ms (c), 4ms (d), 5ms (e), 6ms (f), 7ms 
(g) and 8 ms (h). Co-ordinate axes are inset and the horizon- 
tal bar indicates the scale. 

IV. REFLECTION PROBABILITY OF THE BEC 



In this section, we investigate how the reflection prob- 
ability of the BEC varies with v x and depends on the 
presence or absence of an initial vortex. Defining the 
reflection probability of the condensate as 



R(v x = Axuj x ) 



ly Sx<o I ^ ^ 1 = 71 I 2 dxd y 

J y f x I Vfay,t) | 2 dxdy 



(7) 

we use our numerical solution of the Gross-Pitaevskii 
equation to quantify the reflection probabilities of BEC 
A [R(v x ), dotted curve in Fig. 9(a)] and its irrotational 
equivalent [Ri(v x ), dashed curve in Fig. 9(a)]. Figure 
9(a) shows that as v x decreases, both R(v x ) and Ri(v x ) 
increase. For v x < 7.3mm s _1 , R(v x ) < Ri(v x ) and 
conversely for v x <^ 7.3mm s _1 , R(v x ) > Ri(v x ). This 
crossover is revealed more clearly by the dashed curve 
in Fig. 9(b), which shows r(v x ) = R(v x ) — Ri(v x ). For 
v x ^ 7.3mm s _1 , r(v x ) < 0, but at higher v x , r(y x ) > 0. 

In Fig. 9(a) we compare R(v x ) and Ri(v x ) to 
the equivalent plane wave reflection probability Rp(v x ) 
(solid curve). The difference between Rp(v x ) and Ri(v x ) 
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FIG. 9: (a) Dotted curve: reflection probability, R(v x ), for 
BEC A. Dashed curve: reflection probability, Ri(v x ), for 
BEC Ai. Solid curve: reflection probability, Rp(v x ), for a 
plane wave, (b) Dashed curve: r(v x ) obtained from numeri- 
cal simulations. Solid curve: the expression for r(v x ) obtained 
from Eq. (18) with I 0.007 mm 2 s~ 2 and J 0.19 mm 2 
s -2 . (c) Solid curve: the local velocity distribution function, 
fi(vx — v x , y = 0) for BEC Ai. Dashed curve: the local veloc- 
ity distribution function, f(v x —v x ,y = 0) for BEC A. Dotted 
curve: g(v x - v x , y = 0). 



is due to the finite width of the BEC along the in- 
direction, which means that the atoms start from var- 
ious positions in the harmonic trap and thus reach the 
barrier with a range of incident velocities rather than 
the single velocity of an incident plane wave. We can 
model the spread of velocities by defining the probabil- 
ity Fi(v x -v x )dv x = J y fi(v x -v x , y)dydv x that an atom 
arrives at the barrier with an incident velocity between 
v x and v x + dv x . Here fi(v x - v x ,y)dv x = \Sb((v x - 
v x )/uj x ,y : t = 0)\ 2 dv x / j Vx \V((v x -v x )/uj x ,y,t = 0)\ 2 dv x 
is the local probability that an atom starting from posi- 
tion y arrives at the barrier with a normal velocity com- 
ponent between v x and v x + dv x . For each y, this local 
probability is determined from the initial shape of the 
BEC wavefunction, which defines the initial spatial dis- 
tribution of the atoms and hence the resulting velocity 



distribution as the atoms arrive at the barrier. For each 
y, fi(v x — v x ,y) is symmetrical about v x — v x = owing 
to the symmetry of the initial wavefunction. This can 
be seen from the solid curve in Fig. 9(c), which shows 
fi(v x — v x ,y = 0) for BEC Ai. The reflection probability 
of the BEC can be expressed in terms of its velocity dis- 
tribution and the reflection probability of a plane wave 

by 



POO 

Ri(v x ) = / Fi(v x -v x )R P (v x )dv x . 
Jo 



(8) 



Since Fi(v x — v x ) is sharply peaked at v x = v x we can 
use the Taylor expansion 



Rp(v x ) ~ Rp(v x 



dR P (v x ) 



dv x 



d 2 R P (v x ) 



(v x - v x y 



2dv 2 x 

to write Eq. (8) in the form 

f>00 

Ri(v x ) w Rp(v x ) / Fi(v x -v x )dv. 
Jo 

dR P (v x ) 



(9) 



+ 



dv x 
d 2 R P (v x ) 



Fi(v x - v x )(v x - v x )dv x 



2dvl 



f>00 

I Fi(v x - v x )(v x - v x ) 2 dv x . 

v x JO 

(10) 



Since Fi(v x — v x ) is symmetrical about and its width 
is less than v x for all the v x values considered in Figs. 
9(a,b), Eq. (10) reduces to 



Ri(Vx) « Rp(Vx) + " —g^ 2 



where 



poo 

/ Fi(v x -v x )(v x -v x ) 2 
Jo 



dv x > 0. 



(ii) 



(12) 



It follows from Eqs. (11) and (12) that Ri{v x ) > Rp(v x ) 
if 



d 2 R P (v x 



dvl 

and Ri(v x ) < Rp(v x ) if 

d 2 R P (v x ) 



> 



< 0. 



(13) 



(14) 



Evaluating the integral in Eq. (12) numerically, we find 
that for BEC A u I = 0.007 mm 2 s~ 2 . 

From the shape of the solid curve in Fig. 9(a), we see 

that for v x ^ 7.3 mm s _1 , - _ < 0, whilst for 
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v x <; 7.3 mm s x , 



> 0. Hence, from Eq. 



(11), we expect Ri(v x ) > Rp(v x ), for v x ^ 7.3mm s _1 , 
and Ri(v x ) < Rp(vx), f° r v x ^ 7.3 mm s _1 , which is 
confirmed by comparing the dashed [Ri(v x )] and solid 
[Rp(v x )] curves in Fig. 9(a). 

We now consider a BEC that initially contains a single 
vortex. In this case, Eq. (8) can be generalized to 



R(v x ) 



f 

Jo 



F(v x - v x )R P (v x )dv x 



(15) 



where the velocity distribution for the incident atoms in 
the rotating BEC, F(v x — v x ), equals the velocity distri- 
bution function for an irrotational BEC plus a correction 
G(v x — v x ) that is 



F(v x - v x ) = Fi(v x - v x ) + G(v x - v x ), 



(16) 



where F(v x - v x ) = f y f(v x - v x , y)dy and G(v x - v x ) = 
J d( v x — v Xl y)dy are defined in terms of the local velocity 
distribution functions f(v x —v x ,y) and g(v x —v x ,y) = 
f(v x -v x , y)-fi(v x -v x , y), which are shown in Fig. 9(c), 
for y = 0, by the dashed and dotted curves respectively. 
Since g(v x — v x ,y) and G(v x —v x ) are both symmetrical 
about v x — v x = 0, it follows from an analysis similar to 
that presented in Eqs. (8-10) that 



R(v x ) 



+ 



Ri(v x ) 
d 2 R P (v x ) 



Idvl 



f 



G(v x - v x )(v x - v x ) 2 dv x 



(17) 



which can be rewritten as 
R(v x ) « Ri(v x ) + 

~ Rp(v x ) + 



J d 2 R P (v x ) 



2 dv 2 x 
(/+ J) d 2 R P (v x ) 



dvl 



where 



poo 

J = I G(v x -v x )(v x -v x ) 2 
Jo 



dv x . 



(18) 



(19) 



Unlike /, the sign of J can be either negative or posi- 
tive, depending on the form of G(v x —v x ). To highlight 
this, we initially consider a simple analysis based on the 
form of the local distribution function g(v x — v x ,y — 0) 
[Fig. 9(c)]. For v x w v x (i.e. in the region of the vor- 
tex core) f(v x — v x ,y = 0) [dashed curve in Fig. 9(c)] 
is close to zero, because the atom density falls to zero at 
the vortex core, whereas fi(v x —v x ,y = 0) [solid curve 
in Fig. (c)] is maximal. Consequently, g(v x — v x ,y = 0) 
[dotted curve in Fig. 9(c)] is negative for v x ~v x . Since 
the atoms have moved away from the vortex core, to- 
wards the edges of the BEC, we expect that away from 
the vortex core (where \v x — v x \ ^> 0) f(v x — v x ,y = 0)> 
fi(v x ~v x ,y = 0) and hence g(v x -v x ,y = 0)> 0, as 



can be seen in Fig. 9(c) (dotted curve). Consequently, 
the contribution to J is likely to be positive because the 
integrand in Eq. (19) is largest when (v x —v x ) 2 is large. 

However, if y ^ 0, v x is also perturbed by the circu- 
lation of the BEC around the vortex core. For y > 0, 
v x is increased and for y < 0, v x is decreased. This ef- 
fect further increases the spread of incident velocities at 
the barrier and, due to the (v x — v x ) 2 term in Eq. (19), 
makes J more positive. From our full numerical analysis 
we find that the dominant cause of the shift in the re- 
flection probabilities is the circulation of atoms around 
the vortex. Evaluating Eq. (19) numerically, we find for 
BEC A that J = 0.19 mm 2 s~ 2 > I. 

Since J and / are both positive for the BEC con- 
sidered here, from Eqs. (11) and (18), we expect that 
R(v x ) — Rp(vx) will have the same sign as, but a larger 
magnitude than, Ri(v x ) — Rp(v x ). This is confirmed by 
the curves shown in Fig. 9(a) and highlights the fact 
that the presence of the vortex increases the distribution 
of incident velocities upon the scattering potential. 

In Fig. 9(b), we compare the expression for r(v x ) w 

J d 2 R P (v x ) 
2 dvl 



, obtained from Eq. (18) (solid curve), with 

c 

the numerical values of r(v x ) (dashed curve) obtained 
from our solutions of the Gross-Pitaevskii equation. The 
expression for r(v x ) (solid curve) and the numerical sim- 
ulations (dashed curve) are in good agreement because 
the spread of incident velocities is narrow and hence the 
Taylor expansions used in the derivation of Eqs. (11) and 
(18) are reasonably accurate. For BEC B we find that 
the discrepancy between the numerical r(v x ) values and 
those obtained from Eq. (18) becomes larger. This dis- 
crepancy is not only due to an increase in the spread of 
incident velocities, which introduces higher order terms 
into the Taylor expansion, but also arises from the dis- 
ruption of the BEC upon scattering, an effect which is not 
described in our approximate analysis. Equations. (11) 
and (18) are therefore only valid when no interferential 
or rotational disruption is present in the BEC. 



V. CONCLUSIONS 

We have investigated how BECs with different atom 
densities scatter off a Gaussian potential when the initial 
state is either irrotational or contains a single vortex. 
We find three distinct regimes for the formation of dy- 
namical excitations and hence rotational or interferential 
disruption: (i) at low densities there is no fragmentation 
of the reflected or transmitted components of the BEC, 
irrespective of whether or not there is a vortex in the ini- 
tial state; (ii) at medium densities for which t r /t sr ^ 1, 
rotational disruption occurs in the reflected component 
of a BEC with an initial vortex, but no disruption is ob- 
served when the BEC is initially irrotational; (iii) at high 
densities, there is strong interferential disruption in the 
reflected atom cloud if t c ^ t Sl both in the presence and 
absence of an initial vortex. 
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By considering the velocity distribution of the inci- 
dent atoms, we have derived expressions for the reflection 
probabilities of rotating and irrotational BECs in terms 
of the reflection probability of a single plane wave inci- 
dent on the scattering potential. This analytic approach 
agrees well with our numerical calculations of reflection 
probabilities for BECs scattering from a Gaussian bar- 
rier. It shows that the velocity spread of an irrotational 
BEC causes a positive or negative deviation from Rp(v x ), 
depending on the curvature of Rp(v x ). When a vortex is 
introduced, the circulation further increases the spread 



of incident velocities, leading to an even larger deviation 
from Rp(v x ). 

Finally, we note that with current techniques it should 
be possible to perform experimental tests of our theoret- 
ical predictions, relating to the existence of three distinct 
dynamical regimes and to the effect of a vortex on the 
reflection probability of a BEC impinging on a potential 
barrier. 

This work was supported by the ARC, the EPSRC, the 
Royal Society (London) and the University of Melbourne. 



[1] O. Morsch, J.H. Miiller, M. Cristiani and E. Arimondo, 

Phys. Rev. Lett. 87, 140402 (2001). 
[2] A. S. Arnold, C. MacCormick, and M. G. Boshier, Phys. 

Rev. A 65, 031601 (2002). 
[3] K. Bongs, S. Burger, G. Birkl, K. Sengstock, W. Ertmer, 

K. Rzazewski, A. Sanpera, and M. Lewenstein, Phys. 

Rev. Lett. 83, 3577 (1999). 
[4] A. Giinther, S. Kraft, M. Kemmler, D. Koelle, R. Kleiner, 

C. Zimmermann and J. Fortagh, Phys. Rev. Lett. 95, 

170405 (2005). 

[5] T.A. Pasquini, Y. Shin, C. Sanner, M. Saba, A. Schi- 

rotzek, D.E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 

93, 223201 (2004); T.A. Pasquini, M. Saba, G.-B. Jo, Y. 

Shin, W. Ketterle, D.E. Pritchard, T.A. Savas and N. 

Mulders, Phys. Rev. Lett. 97, 093201 (2006). 
[6] R.G. Scott, A.M. Martin, T.M. Fromhold and F.W. 

Sheard, Phys. Rev. Lett. 95, 073201 (2005). 
[7] M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 

Hall, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 

83, 2498 (1999). 
[8] K. W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 

ibard, Phys. Rev. Lett. 84, 806 (2000). 
[9] F. Chevy, K. W. Madison, and J. Dalibard, Phys. Rev. 

Lett. 85, 2223 (2000). 
[10] C. Raman, J. R. Abo-Shaeer, J. M. Vogels, K. Xu, and 

W. Ketterle, Phys. Rev. Lett. 87, 210402 (2001). 
[11] R.G. Scott, A.M. Martin and T.M. Fromhold, Phys. Rev. 

A 69, 063607 (2004). 
[12] J.J. Garca-Ripoll, G. Molina- Terriza, V. M. Prez-Garca 

and L. Torner, Phys. Rev. Lett. 87, 140403 (2001). 
[13] I. Josopait, L. Dobrek, L. Santos, A. Sanpera and M. 

Lewenstein, Eur. Phys. J. D. 22, 385 (2003). 
[14] Y. Shin, M. Saba, T.A. Pasquini, W. Ketterle, D.E. 

Pritchard and A.E. Leanhardt, Phys. Rev. Lett. 92, 

050405 (2004). 

[15] M. Greiner, O. Mandel, T.W. Hansen and I. Bloch, Na- 
ture 419, 51 (2002). 



[16] R.G Scott, D.A.W . Hutchinson and C.W. Gardiner, 

cond-mat/0608135 (unpublished). 
[17] M.L. Chiofalo, S. Succi and M.P Tosi, Phys. Rev. E, 62, 

7438 (2000). 

[18] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P 
Flannery, Numerical Recipes, the Art of Scientific Com- 
puting (Cambridge University Press, Cambridge, 1994). 

[19] We calculate the expectation values of the positions of the 
transmitted and reflected atom clouds when they come 
to rest at t = tv/cu x as 

L _ I y L>o \y( x >y, t = */ux) | 2 xdxdy 
T ~ f y f x \*(x,y;t) \idxdy 

and 

L _ Iy L<o I ^ ( x ' = I 2 xdxd y 

J y J x \V(x,y;t) | 2 dxdy 

[20] Since the interactions in BEC C play a dominant role, 
the diameter of the condensate is 50/xm and hence it is 
only possible to consider Ax > 25/rni. For this reason we 
only consider Ax = 30/xm for system C. 

[21] R.G. Scott, A.M. Martin, T.M. Fromhold, S. Bujkiewicz, 
F.W. Sheard and M. Leadbeater, Phys. Rev. Lett. 90, 
110404 (2003). 

[22] R.G. Scott, A.M. Martin, S. Bujkiewicz, T.M. Fromhold, 

N. Malossi, O. Morsch, M. Cristiani and E. Arimondo, 

Phys. Rev. A. 69, 033605 (2004). 
[23] L. Fallani, L. De Sarlo, J.E. Lye, M. Modugno, R. Saers, 

C. Fort and M. Inguscio, Phys. Rev. Lett. 93, 140406 

(2004). 

[24] Z. Dutton, M. Budde, C. Slowe and L.V. Hau, Science 

293, 663 (2001). 
[25] Clearly, for a symmetrical condensate, j dv x F(v x — v x )> 

0. 



